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Abstract 

A model is presented for obtaining the step formation energy for 
metallic islands on (111) surfaces from Monte Carlo simulations. This 
model is applied to homo (Cu/Cu(lll), Ag/Ag(lll)) and heteroepi- 
taxy (Ag/Pt(lll)) systems. The embedded atom method is used to 
represent the interaction between the particles of the system, but any 
other type of potential could be used as well. The formulation can 
also be employed to consider the case of other single crystal surfaces, 
since the higher barriers for atom motion on other surfaces are not a 
hindrance for the simulation scheme proposed. 

Keywords: Step formation energy, Monte Carlo simulation, sub- 
monolayers. 



1 Introduction 

The formation of metal islands on a metal surface represents an important 
stage in the growth of a phase in the case of homoepitaxy, or the appearance 
of a new one in the case of heteroepitaxy. In any case, this topic has interested 
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whole generations of theoreticians and experimentalists, and is currently the 
subject of intense research. Just trying to cite relevant work in this field 
would require by itself a whole extensive review, so we just mention briefly 
the work which may be connected directly to the present letter. 

Concerning homoatomic systems, recent experimental work has allowed 
the accurate determination of the equilibrium shape of two dimensional metal 
islands on different single crystal surfaces The analysis by the inverse 
Wulff construction allowed the determination of the ratio of the free energies 
per step length of the A- and B- type steps on Cu(lll) and Ag(lll) surfaces, 
as well as the angle dependence of the step free energy for Cu(100), Cu(lll) 
and Ag(lll). The analysis of the experimental data using Ising models 
worked quite well for the (111) surfaces but were not helpful in general for 
the (100) surfaces. On the other hand, an alternative method for determining 
the kink energy from Arrhenius plots was developed, which did not rely on a 
specific model for the interactions between the atoms but on thermodynamic 
arguments. An important quantity also emerging from these studies was the 
step energy per atom, which was obtained for straight and 100 % kinked 
steps, thus yielding a detailed picture of the energetics of these systems. 

In the case of heteroatomic systems, the interatomic interactions and the 
adsorbate/substrate misfit determine the structure and the growth mech- 
anism of the monolayer and subsequent adlayers||. Although the experi- 
mental data on these systems is extensive, information concerning the step 
formation energy is scarce. This quantity has been estimated from electro- 
chemical measurements by Xia et al [[| for the system Cu/Au(lll), finding 
a value of 0.4-0.5 eV/(step atom). As we shall see below, in the case of de- 
position of Ag on Pt(lll) [|J a rough estimation of this quantity can also be 
made from the temperature at which the fragmentation of islands occurs. 

It is the purpose of the present work to formulate a model for the estima- 
tion of step formation energies from Monte Carlo simulations, and perform 
calculations for some typical systems. Our model includes the natural relax- 
ation of the system due to temperature effects, which coupled to the misfit 
between substrate and adsorbate may produce some singular effects in cer- 
tain systems. Due to the simplified potential employed, we do not expect to 
make accurate quantitative predictions, but we think the present results will 
motivate further research in the field. The systems studied were Cu/ Cu(lll), 
Ag/Ag(lll) and Ag/Pt(lll). 
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2 The Drop model 



In the present calculations we focus on flat islands that are only one atom 
thick. Both the energy of adsorption per atom in the bulk of the island and 
edge may be of interest. In order to separate the two contributions, we have 
adapted the drop model || that is often used for the calculation of surface 
energies to the two-dimensional situation. This model permits to obtain the 
binding energy of the atoms in the monolayer u s and the binding energy of 
the atoms at the border of the island u b from the adsorption energy per atom 
u t , the number of border atoms N b and the number of inner atoms N s that 
can be obtained for the Monte Carlo simulation. 

Let Ut be the excess of energy of a system composed of an island adsorbed 
on a substrate surface M: 

Ut = (Ua/m — U M ) 

where Ua/m is the energy of the adsorbed-substrate system and Um is the 
energy of the substrate without adsorbate. 

If we consider the adsorbed submonolayer like a 2-D drop, Ut can be 
expressed in terms of the edge atoms N b and the number of atoms which 
make up the surface N s as: 

U T = u s N s + u b N b (1) 
Let iV be the total number of atoms in the island: 

N = N S + N b 

If we replace N s from this equation into equation [I] , we obtain: 

U T = u s N +(u b -u s )N b 
And dividing by iV we get: 

u t = — = u s + (u b -u s )— (2) 

Thus, a linear relationship is expected between the average values < u t > 
and <! ^ > for simulations employing islands of different sizes. The u s and u b 
were obtained using least squares to calculate the ordinate and slope of the 
straight line respectively. 
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The ordinate u s , corresponding to the limit — ► should approach the 
binding energy of a monolayer. The slope (it& — u s ) is the border formation 
energy per atom. This should always be a positive quantity and its amount 
indicates the energy change required to bring an atom from the surface to 
the edge of the island. 



3 Model and Computations 

The substrate surface was assumed to be smooth without considering steps 
or reconstruction. The substrate employed in Monte Carlo simulations was 
a 4-layer slab with 480 atoms per layer to represent (111) crystalline surface. 
The atoms in the first and the second planes of the substrate located within a 
circular area containing ca. 250 atoms were allowed to vary the positions dur- 
ing the simulation. The remaining substrate atoms were fixed to their bulk 
equilibrium configuration to emulate the presence of a semi-infinite crystal. 
On this substrate, submonolayers of an initially circular shape containing 
between 12-150 atoms were adsorbed to emulate islands of different sizes. 

The potentials employed were given by the embedded atom method (EAM) 
which takes into account many body interactions characteristic of the metallic 
systems |J,0. The total energy is expressed in terms of two contributions: 

N x 

i=l i j^i 

where p^i is the host electronic density at i atom position, is the 

embedding energy and fajfaj) is the core-core pair repulsion between i and 
j atoms separated at a r^- distance. In all cases we used the parametrization 
of Foiles et al 0. 

All the Monte Carlo simulations were performed in the (NVT) ensemble 
at 300K. We considered 5000 equilibration steps followed by 15000 production 
steps. 

In order to let the adatoms overcome high energy barriers, like those for 
the displacement on the surface or the detachment from a cluster, we allowed 
these adatoms to perform long jumps, with displacements Ar: 

Ar = Hi's! + n 2 's2 
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where n\ and n 2 are integers and s% and s 2 are the primitive vectors of a 
two-dimensional Bravais lattice. These "long" jumps are important to get a 
proper equilibration of the system. "Short" jumps, as usually employed in 
Monte Carlo simulations where the positions of the particles are varied con- 
tinuously, are also present in our studies to describe the vibrational motion 
of the adsorbate atoms in the neighborhood of their equilibrium positions. 
This simulation method was employed in a previous work || to study Pd 
submonolayer growth on different monocrystalline surfaces, yielding results 
in reasonable agreement with experimental data of literature. 



4 Results and Discussion 

The equilibrium shapes of two dimensional islands have been the subject of 
extensive research in both the experimental and the theoretical fields. As 
mentioned above, a recent work of Giesen et al has analyzed this problem in 
detail [0. While our current computer capabilities do not allow the perfor- 
mance of simulations able to predict the shapes of island with sizes in the 
order of tenths of nanometers, we shall tackle the problem of predicting step 
formation energies. 

We performed simulations for the homoatomic systems Cu/Cu(lll) and 
Ag/Ag(lll) and for the heteroatomic Ag/Pt(lll) system for a series of is- 
lands of different sizes. The values of < u t > and <I ^ > were obtained from 
the simulation production steps for the different submonolayer sets and are 
plotted in Figure 1. We also show there the < u t > values obtained from 
a simulation with a monolayer. We have plotted these monolayer energies 
assigning them the value <I ^ > = 0. For the larger N employed in our simu- 
lations, we observe a linear behavior for the three systems studied. A least 
square fitting in the range 0.2 < <J ^ b> < 0.4 using eqn. (H) delivered the 
u s and Ut> reported in Table 1. This set corresponds to the larger islands 
simulated. For Cu/Cu(lll) and Ag/Ag(lll) the step formation energies 
Ub — u s are in the order of 0.2-0.3 eV, which are close to the values found by 
Giesen et al. 0. In the case of the Ag/Pt(lll) system, an unusually low 
step formation energy is found. This system presents an important positive 
misfit (4.3%) which generates a strong compressive stress in the pseudomor- 
phic (lxl) monolayer ||. Roder et al||] have shown that when islands are 
grown at low temperatures and subsequently annealed, a remarkable mor- 
phology transition is found around 300-350 K. At these temperatures, island 
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fragmentation occurs via edge roughening and kink evaporation, so that at 
400 K all large Ag islands are eventually disintegrated. Generally speaking, 
island disintegration should be expected to take place when the step forma- 
tion energy is of the order of kT. Our estimation of step formation energy 
of ca. 0.085 eV, small compared with those of other systems, delivers a dis- 
integration temperature (990 °K) which results far above the experimental 
value. While this disagreement could be roughly attributed to the approx- 
imate nature of the potential employed, close observation of figure 1 shows 
that the linear extrapolation of the < u t > vs <r ^ > curve to <! ^ > — > shows 
the strongest deviation with respect to the < u t > of the monolayer for the 
Ag/Pt(lll) system. Thus, the question arises if some of the approximations 
involved in the present model are responding bad for this system. In this 
sense, the extrapolation to <I ^ > — > to obtain the ordinate u s from equa- 
tion (0) implies the assumption that the island is roughly composed of two 
sorts of atoms, say "border " and "inner " atoms. In order to learn how fast 
the atoms inside the island approach the behavior of atoms in a monolayer 
( "inner " atoms), we have plotted in figure 2 the distance between nearest 
neighbors at the center of the monolayer as a function of the total number of 
atoms in the island. In the case of the system Ag/Ag(lll), we see that for 
relatively small islands the atoms at the center already reflect the behavior 
of the monolayer. In the case of the Cu/Cu(lll) system, the approach to 
the monolayer behavior is more sluggish, but is practically reached for the 
large island in our simulation. On the other hand, we see that in the case of 
the Ag/Pt(lll) system the behavior of the inner atoms is far from that of 
Ag atoms in a (lxl) monolayer on Pt(lll). These results indicate that small 
Ag islands on Pt(lll) are considerably expanded with respect to the (lxl) 
monolayer, becoming more and more compressed with growing island size, 
probably reaching finally a (lxl) structure with the substrate for sizes that 
are beyond our simulation conditions. In this way, the binding energy of the 
atoms in the monolayer is not accurately estimated from the extrapolation 
in figure 1. On the other hand, the w s values can be properly obtained from 
a simulation employing a (lxl) adsorbed monolayer, and the linear fitting of 
the points in figure 1 can be made by involving Ub as unique free parameter. 
Table 2 shows u s values stemming from monolayer simulations, as well as the 
Ub obtained from linear fittings taking u s from these monolayer simulations. 
The corresponding step formation energies Ub — u s are also reported there. 
Comparison with Table 1 indicate that although the extrapolated ^values 
were reasonable, important changes occur in the step formation energy, es- 
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pecially in the case of Ag/Pt(lll), where w fe — u s has strongly decreased. 

A last point that we have checked is the influence of relaxation in the 
calculated adsorption and step formation energies. In Table 3 we show the 
quantities obtained from static calculations, where the island and surface 
atoms are not allowed to relax. As in the previous case, the most drastic 
differences are found in the case of the step formation energy of the system 
Ag/Pt(lll). 

For all the island sizes of Ag/Pt(lll) considered in this work, we found 
an outwards relaxation of the edge atoms of the islands, which confirms the 
statement of Bromann et al |10[ concerning the fact that compressive stress 
can be relieved by this way. However, we did not find pseudomorphic growth 
in our islands, but rather expanded structures which become progressively 
compressed as the islands increase their size. Simulations using periodic 
boundary conditions showed that the (lxl) structure was stable under our 
simulation conditions. Thus, the picture that emerges from the present cal- 
culations is that the increased size of the islands leads to progressively com- 
pressed structures of Ag on Pt(lll) until they reach the (lxl) structure of 
the monolayer. The actual behavior of the experimental system is richer, 
with the appearance of partial dislocations, yielding domain walls between 
fee and hep stacking |l(J. In this respect, simulations with a much larger 
number of atoms are desirable in order to investigate if the EAM potential 
is able to reproduce such subtle features. 

In contrast to the Ag/Pt(lll) system, the other systems present inwards 
relaxation at the edge, following the natural expectation from bond-order 
conservation analysis, atoms at the edge are less coordinated, and they in- 
crease the strength of the binding with each other. 

We can summarize the present results by stating that the model we have 
presented permits to obtain the step formation energy for metallic islands on 
(111) surfaces in a straightforward way for homo and heteroepitaxy systems 
from Monte Carlo simulations. We have used the embedded atom method to 
represent the interaction between the particles of the system, but any other 
type of potential could be used as well. The model can also be used for other 
single crystal surfaces, since the higher barriers for atom motion on other 
surfaces are not a hindrance for the simulation scheme proposed. 
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5 Tables 



System 


u s [eV] 


u b [eV] 


u b - u s [eV] 


Cu/Cu(lll) 


-3.509±0.001 


-3.194±0.004 


0.315±0.003 


Ag/Ag(lll) 


-2.817±0.001 


-2.583±0.004 


0.234±0.003 


Ag/Pt(lll) 


-3.167±0.002 


-2.997±0.008 


0.085±0.006 



Table 1: Binding energy per atoms in the monolayer u s , binding energy 
per atom at the edge of the island u b and step formation energy per atom 
Ub — u s for different systems. These values were obtained from linear fitting 
of the points in figure 1 in the range 0.2 < <J Y r < 0.4. The parameters fit- 
ted were u s and u\> — u s according to equation ([J). The Monte Carlo/EAM 
simulations were performed at 300 K. 



System 


u s [eV] 


u b [cV] 


u h - u s [eV] 


Cu/Cu(lll) 


-3.498 


-3.216±0.002 


0.283±0.002 


Ag/Ag(lll) 


-2.808 


-2.600±0.001 


0.208±0.001 


Ag/Pt(lll) 


-3.150 


-3.118±0.002 


0.032±0.002 



Table 2: Binding energy per atom in the monolayer u s , binding energy per 
atom at the edge of the island w^and step formation energy per atom Ub — u s 
for different systems. The u s were obtained from monolayer simulations. The 
Ub values were obtained from linear fitting of the points in figure 1 in the 
range 0.2 < < 0.4 according to equation (g). The Monte Carlo/EAM 

simulations were performed at 300 K. 
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System 


u s [eV] 


u b [eV] 


u b - u s [eV] 


Cu/Cu(lll) 


-3.555±0.003 


-3.21±0.01 


0.350±0.007 


Ag/Ag(lll) 


-2.861±0.002 


-2.597±0.008 


0.264±0.006 


Ag/Pt(lll) 


-3.164±0.002 


-3.001±0.003 


0.163±0.001 



Table 3: Binding energy per atom in the monolayer u s , binding energy 
per atom at the edge of the island u b and step formation energy per atom 
u b —u s for different systems, as obtained from static EAM calculations. These 
values arise from a linear fitting of the adsorption energy per atom in the 
range 0.2 < < ^ f b> < 0.4. The parameters fitted were u s and u b — u s according 
to equation (Q). 
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Figure 1: Adsorption energy per atom < u t > vs fraction of border 
curve for different adsorbate(island)/substrate(lll) systems. 
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Figure 2: Nearest neighbor distance at the center of the island as a function 
of the total number of island atoms for the systems considered in this work. 
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